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0\ ' Abstract 

o : 
\o . 

■ This paper reviews the apphcation of a novel methodology for analysing the 

. isotropy of the universe by probing the alignment of local structures in the CMB. 

Qh| The strength of the proposed methodology relies on the steerable wavelet filtering of 

O ' the CMB signal. One the one hand, the filter steerability renders the computation of 

-fi . the local orientation of the CMB features affordable in terms of computation time. 

^ I On the other hand, the scale-space nature of the wavelet filtering allows to explore 

^ ' the alignment of the local structures at different scales, probing possible different 

phenomena. We present the WMAP first-year data analysis recently performed by 
the same authors (Wiaux et al.), where an extremely significant anisotropy was 
I found. In particular, a preferred plane was detected, having a normal direction 

with a northern end position at {0,ip) = (34°, 331°), close to the northern end of 
the CMB dipole axis. In addition, a most preferred direction was found in that 
plane, with a northern end direction at {9,ip) = (71°, 91°), very close to the north 
ecliptic pole. This result synthesised for the first time previously reported anomalies 
identified in the direction of the dipole and the ecliptic poles axes. In a forthcoming 
paper (Vielva et al.), we have extended our analysis to the study of individual 
frequency maps finding first indications for discarding foregrounds as the origin of 
the anomaly. We have also tested that the preferred orientations are defined by 
structures homogeneously distributed in the sky, rather than from localised regions. 
We have also analysed the WMAP 3-year data, finding the same anomaly pattern, 
although at a slightly lower significance level. 
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1 Introduction 



The release of the Wilkinson Microwave Anisotropy Probe (WMAP) first- 
year data in 2003 (U and the WMAP 3-year data in 2006 Q) have implied a 
huge activity in the Cosmic Microwave Background (CMB) community during 
the last years. Besides the impressive accuracy in the determination of the 
cosmological parameters given by the WMAP data analysis (0; jij) as well 
as by combining other cosmological data sets (e.g. (j^)), a large effort was 
concentrated on the very challenging issue of probing the stationary Gaussian 
random nature of the CMB predicted by the standard inflationary Big-Bang 
model. This is a capital topic, since the determination of the cosmological 
parameters that define the concordance cosmological model relies on those 
assumptions. 



Strong evidences (at confidence levels above 99%) supporting the non-stationarity, 
or statistical anisotropy, of the CMB signal were established. Analyses based 
on N-point correlation functions JE local curvature (0), local power spec- 
tra ^ 11), and bispectra suggest a north-south asymmetry max- 
imised in a coordinate system with the north pole at {6, ip) = (80°, 57°) in 
Galactic co-latitude 9 and longitude ip, close to the north ecliptic pole lying 
at {9,ip) = (60°, 96°). Analyses of multipole vectors, angular momentum dis- 
persion, as well as azimuthal phases correlations find an anomalous alignment 
between the low / multipoles of the CMB, suggesting a preferred direction 
around (6*, ip) = (30°, 260°), near the ecli ptic plane and close to the axis of the 
dipole lying at {6, ip) = (42°, 264°) (Qlii [li H Hlil El B- Galactic 
north-south asymmetries are also found in the analysis of the kurtosis and 
the area of the wavelet coefficients of the CMB data (j2l|). These are mainly 
due to a very cold spot (the so-called cold spot) in the southern hemisphere 
(H H H H Hi) and Q in this volume). First results with the angu- 
lar pair separation method, which probes the statistical isotropy both in real 
and multipole space, also seem to support those results (|27ll. In opposition 
to the previous results, bipolar power spectra analyses are consistent with no 



311). Finally, 



violation of the statistical isotropy of the universe ()28l : l29l : |3C 
theoretical models for an anisotropic universe are being studied to account for 
the observed effects Q H S S H . 



The method we introduced in ()38l ) probes the statistical isotropy of the CMB in 
a novel way: the analysis of the alignment of structures of the CMB signal. By 
searching for the directions towards which the local CMB structures are mostly 
oriented, preferred directions in the sky may be defined. The method is based 
on the steerable wavelet decomposition of the CMB. On the one hand, the filter 
steerability renders affordable in terms of computation time the calculation of 
the local direction of the CMB structures at each position on the sky. On the 
other hand, the scale-space nature of the wavelet filtering allows the probe of 
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An even signal is abtaiaeii 



Fig. 1. Graphical illustration of the proposed methodology. The positions on the 
sphere towards which the local CMB structures are aligned are defined as those lying 
on the great circles that pass by the structures and are parallel to their orientations. 
The signal defined in this way is even (in Cartesian coordinates) by construction. 
This construction probes preferred directions in the sky towards which the CMB 
features are mostly oriented. As such, it represents a powerful analysis of the hy- 
pothesis of the universe isotropy. This figure is a modification of one from the Max 
Tegmark web site (|http://space.mit.edu/home/tegmark/index.html|. 

the local CMB structures alignment as a function of the scale considered. Our 
analysis synthesised for the first time previously reported anomalies related 
with the dipole and the ecliptic poles axes. In this contribution we review our 
first work (jssh as well as some new results from our forthcoming paper (E^, 
relative of the detected anomaly. In Section 2 the details of the methodology 
are given. The application to WMAP data is detailed in Section 3. Finally, 
the Conclusions are in Section 4. 



2 Methodology 

Under the cosmological principle assumption, and in the framework of the 
basic inflationary models, the CMB can be interpreted as a stationary Gaus- 
sian random field on the sphere. The hypothesis of stationarity, or statistical 
isotropy, implies that no preferred orientation is expected. We propose to 
probe the statistical isotropy of the CMB by studying the orientation of its 
structures. The number of times a given direction in the universe is seen by the 
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Fig. 2. Mollweide projection of the first member ^'02 ^ of the second Gaussian deriva- 
tive (2GD) wavelet basis at position loq = (6*0, v?o) = (7r/2,0) and scale a = 0.19. 

local CMB features provides a unique way to determine whether there exists 
or not any preferred direction on the sky towards which the CMB structures 
are unexpectedly aUgned. 

An illustration of the proposed methodology is given in Figure 1. Let us focus 
on the two particular CMB structures in the middle of the two dark circles. All 
the positions on the sphere lying on the great circles that pass by the selected 
structures and are parallel to their orientations, are said to be seen by the 
corresponding structure. As it is indicated in the Figure, there is a particular 
position in the sky that is seen by both structures. Once this procedure is done 
for all the CMB structures, we end up with a signal on the sphere, allowing for a 
large number of statistical analyses (e.g. angular power spectrum, geometrical 
statistics, Minkowski functionals, ...). This signal is, by construction, even 
(in Cartesian coordinates), since it is obtained by analysing the direction of 
the features, without any notion about the specific sense. 

The method is based on the steerable wavelet decomposition of the CMB 
signal. The reason for working in wavelet space is twofold. First, the wavelet 
filtering provides the natural framework for a multi-scale analysis of the local 
CMB structures. This opens the possibility to probe the statistical isotropy 
of the CMB at different angular resolutions, which is crucial since a possible 
source of anisotropy may be due to structures of a particular size. Second, 
steerable filters render the computation of the orientation of the features at 
each position on the sky affordable in terms of the computation time 
Steerable filters were firstly introduced on the plane by (|4^ and they have 
been recently extended to the sphere by (|43h. A non-axisymmetric filter is 
steerable if any rotation by x ^ [0, 27r) around itself may be expressed in 
terms of a finite linear combination of non-rotated basis filters. The reader 
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is kindly referred to the previous references for a detailed description of the 
corresponding properties. The second Gaussian derivative (2GD) is a steerable 
wavelet filter, defined by the three following basis filters, at any scale a > 0: 



^9|,a {0, if) 



y 




g-2tan2(6»/2)/a2 



g-2tan2(6»/2)/a2 



g-2tan2(e/2)/a2 



sin {2ip) tan^ -, 





4 ^ 6 



where the position u at each point on the sphere is identified by its co-latitude 
6 e [0, vr] and longitude if G [0, 27r). In Figure 2 we plot the first member of the 
basis, a, centred at position uq = (9q, ipo) = (7r/2, 0), and scale a = 0.19. 



3 Application to WMAP data 

We have applied the methodology outlined in Section 2 to the WMAP data. 
The results, obtained on the WMAP first-year data, where published in L3 8.) . 
In a forthcoming paper (jsol ). we extend the first analysis and present the 
results obtained for the three-year data. In this Section we summarise the 
most important results. 

3.1 Data and simulations 

Since we are interested in probing the CMB statistical isotropy, we need to 
deal with a CMB map as free of foreground emission and unknown systematics 
as possible. Many different CMB maps (all them obtained from the WMAP 
data) have been proposed and used during the last years. We considered the 
map defined by ^44") in the WMAP Gaussianity paper. This map (the so-called 
template fits map) is generated as a noise-weighted combination of the eight 
temperature maps at the Q, V and W frequency channels, previously cleaned 
from foreground emission by performing a joint multi-template fitting (see 
dish for details). After the combination, the so-called KpO mask is applied to 
account for diffuse Galactic emission close to the Galactic plane as well as 
for discarding bright point sources. Finally, the residual monopole and dipole 
are subtracted outside the mask. The resultant map is presented in Figure 3. 
The WMAP data are given in the HEALPix pixelization (|4^ at resolution 
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Fig. 3. The analysed CMB map obtained from the WMAP first-year data. This map 
is obtained by following the template fits method described in (|45l ) and firstly used 
by |i3) in the WMAP Gaussianity analysis. 

parameters Nsi^e = 512 and Nside = 1024. However, we have performed our 
analysis on maps downgraded at Ngi^e = 32, where the signal-to-noise ratio 
is extremely high. At that resolution, the sphere is defined by Npi^j. = 12288 
pixels, with an angular resolution of 1.8°. 

In order to confront the obtained results with the expected isotropic behaviour, 
we have performed 10000 simulations. For each simulation, a stationary Gaus- 
sian CMB realisation is generated from the theoretical angular power spectrum 
Ce computed with CMBFAST ()47l ). using the cosmological parameters given 
by the WMAP best-fit cosmological model. Each receiver is simulated by con- 
volving the simulated map with the corresponding beam window function and 
by adding an anisotropic Gaussian noise realisation of the particular noise 
level. As for the data, the eight maps are linearly combined using weights 
defined by the noise amplitude per pixel. The KpO mask is applied, before 
removing the residual monopole and dipole. 



3.2 Results 



We have explored the CMB statistical isotropy at twelve different scales. For 
a given scale a of the wavelet, the angular size of the 2GD wavelet is defined 
as twice the dispersion of the corresponding Gaussian. The twelve scales con- 
sidered, correspond to angular sizes ranging from 5 to 50°. By convolving the 
CMB map with the wavelets, the pixels close to the mask are going to be 
very affected by the zero value of the mask. For that reason, at each scale, an 
extended exclusion mask Ma is defined in order to avoid highly contaminated 
pixels close to the mask. For each pixel outside the mask in the sky and for 
each scale, the orientation for which the wavelet coefficient is maximum in 
absolute value is selected and the corresponding absolute value is retained. 
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Fig. 4. Top panel: Map of total weights Da^ (cu) resulting from the WMAP first-year 
data at scale 03 corresponding to an angular size of 8.3°. Bottom panel, the same 
map, but expressed in (Tag units (estimated from the 10000 simulations). 

This procedure selects the local wavelet orientation that best matches the 
orientation of the local structure of the signal at each point. 

As explained in (jiih. the total weight Da{uj) is a weighted measure of the 
number of times a given pixel uj in the sky is seen by local features at a given 
scale a in the original signal. It is defined as: 

Da{u:) = - Y: \WiU\xoU^),a)l (1) 

^ c=l 



where W^{uJo,Xo{^o):(^) is the maximum value of the wavelet coefficient at 
the position ujq (that is achieved at orientation xo ^ [0,27r)) and the factor 
A = LNpj^Y^woiMa 1^* (i^O! Xo(i^o)) of normalisation defines a mean total 
weight in each direction equal to unity for isotropic CMB simulations without 
mask: A''"^ X^ojgs -Da('^) = 1- The quantity L = ANgi^e stands for the number 
of points on a great circle on a HEALPix grid. 
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Fig. 5. The 20 anomalous directions (pairs of opposite points) with total weights 
higher than in any of the 10000 simulations are shown. They define a preferred 
plane with a perpendicular northern end position at {9,ip) = (34°, 331°), very close 
to the CMB dipole axis at {0,ip) = (42°, 264°) and to the so-called axis of evil. 
A preferred direction in the preferred plane is identified with a northern end po- 
sition at {0,ip) = (71°, 91°), which is extremely close to the north ecliptic pole at 
(9,^) = (60°, 96°). 

In the top panel of Figure 4 we present the signal Da,_^ (uj) obtained from the 
analysis of the WMAP first-year data. The angular size corresponding to the 
scale as is 8.3°. It is clear that the zone around the Galactic plane is less 
seen than the rest of the sky. It is an obvious result, since a pixel has more 
probability to be seen from nearby pixels than from distant ones. Hence the 
exclusion of pixels by the mask implies that, even in an isotropic universe, 
the pixels close to the mask are less seen than the others. Obviously, what 
is much more interesting is to represent the previous signal as normalised, 
at each point, by the standard deviation (sigma) estimated from the 10000 
simulations, i.e. in (Tq units. This normalisation notably cancel out all the 
spurious anisotropies, like the one due to the mask. This map is presented in 
the bottom panel. As it can be seen, there is a very particular pattern with 
many great circles crossing the equatorial poles. Obviously, the distribution 
of Das i^) Gaussian, and, therefore, the number of sigmas is not a 

direct measurement of the probability of each direction, although it is quite 
useful for a visual illustration. In any case, the probability associated with 
each direction in the sky (pairs of opposite points) can be easily computed 
from the simulations. In particular, there are 27 directions that present values 
of {u) anomalous at 99.99% (i.e., at maximum, only one simulation has a 
larger value). Even more, just one simulation has at least 27 directions above 
99.99% at scale ag. 

Among the previous 27 directions, 20 present values of {u) larger than in 
any of the 10000 simulations. These directions (Figure 5) are also identified 
as being anomalous at more than 99.99% (a 100% statement being prohib- 
ited by the finite size of our sample of simulations). These anomalous direc- 
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Fig. 6. The pixels that are ahgned towards any of the 20 directions anomalous at 
more than 99.99% in the WMAP first-year data are shown (in dark red). They are 
homogeneously distributed in the sky, which seems to indicate that the source of 
the detected anomaly is not localised in a specific region of the sky. 

tions are shown in Figure 5. They are mostly concentrated in two clusters of 
4 and 16 directions. The mean direction^ of the latter has a northern end 
at {9, if) = (71°, 91°), which is extremely close to the north ecliptic pole at 
{9,ip) = (60°, 96°). But even more, the two clusters seem to lie on a great 
circle, hence defining a preferred plane in the sky. The direction perpendicu- 
lar to this preferred plane has a northern end position at {6, if) = (34°, 331°) 
which is very close to the CMB dipole axis at [6, ip) = (42°, 264°) and to the 
so-called axis of evil. Let us emphasise that also at scale 04 (around 10° of 
angular size) the CMB presents the same anomaly pattern. In this case there 
are 11 anomalous directions with a total weight Da^ (u) larger than the one 
obtained with any of the 10000 simulations. 

In a forthcoming paper (jl^, we investigate possible origins for this anomaly. 
For instance, one interesting question is to investigate whether this anomalous 
detection is due to some structures spatially localised or, on the contrary, the 
origin for such anomaly is homogeneously distributed in the sky. To answer 
that question, we have plotted in Figure 6 the pixels that are aligned towards 
any of the 20 anomalous direction. One can infer that the latter hypothesis is 
the most plausible one. 

We are also analysing whether the origin of the anomaly could be caused 
by a residual foreground contribution on the CMB map. If so, one should 
expect to find a frequency dependence for the anomaly, as the foregrounds 
themselves are frequency dependent. In Figure 7 we plot the Da.^ [u) signal, in 
(Ta3 units, obtained for each of the Q, V and W WMAP frequency channels. 
The pattern already shown in Figure 5 is as well observed on each one of the 
analysed frequencies, indicating that the anomaly seems to be not caused by 

^ The mean direction is obtained by weighting each direction by its corresponding 
total weight. 
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Fig. 7. The total weight Da^ (lo) are plotted in units for the Q (top), V (middle) 
and W (bottom) WMAP frequency maps. The obtained pattern is the same for all 
of them and is the same as the one already shown in Figure 5. This fact seems to 
discard the foregrounds as responsible for the anomaly. 

foregrounds. 

Finally, the preliminary analysis of the WMAP 3-year data also seems to 
indicate the same structure of anomalies at the same scales. For comparison 
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Fig. 8. Map of total weights Da-^ (a;) resulting from the WMAP 3-year data at 
scale as in flag units. The same pattern as in the first-year data analysis is observed, 
although the significance of the anomaly has slightly decreased. 

with the resuhs for the first-year data (Figure 4, bottom panel), in Figure 8 we 
present the (a;) in units for the 3-year data. The same pattern can be 
recognised in both data sets, although for the 3-year data the significance of 
the anomaly seems to be slightly smaller than in the first-year data: whereas 
in the first release 27 positions were anomalous at 99.99% (i.e., at maximum, 
only one simulation has a larger value), in the second one, only 15 positions 
had a total weight with the same significant. The number of simulations having 
at least 15 positions at that significance or above is 42. 



4 Conclusions 

This contribution reviews a novel methodology for probing the universe isotropy 
through the statistical isotropy of the CMB. It relies on the analysis of the 
alignment of local features in the CMB proposed by the authors in (j38h . The 
power of the method is based on the use of steerable wavelets. On the one 
hand, steerability makes feasible the computation of the local orientation of 
the CMB features which is otherwise unaffordable in terms of computation 
times. On the other hand, the steerability renders accessible the computation 
of the local orientation of the CMB features, which is otherwise unaffordable 
in terms of computation times. On the other hand, the scale-space nature of 
the wavelet scheme allows to explore the alignment of the local structures at 
different scales, probing possible different phenomena. 

In the first application of the method (|38h to the WMAP first-year data, 
we found an extremely significant anisotropy. 20 anomalous directions were 
identified in the sky at a very high confidence level (> 99.99%) at a scale 
corresponding to an angular size of 8.3°, corresponding to a multipole range 
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between i = 11 and i = 27. The 20 anomalous directions where located in 
two clusters of 4 and 16 directions respectively. The mean position of the 
biggest one is at {0,ip) = (71°, 91°), very close to the north ecliptic pole 
{6, (p) = (60°, 96°). In addition, the two clusters are aligned in the same plane, 
with a northern end normal direction pointing towards {0,(p) = (34°, 331°), 
very close to the CMB dipole axis {6, (p) = (42°, 264°) and the so-called axis 
of evil. This result synthesised for the first time previously reported anomalies 
identified towards the dipole and the ecliptic poles axes. 



We are currently analysing the origin of such an anomaly (j39r i as well as 
its significance in the WMAP 3-year data. The second WMAP data release 
provides a pattern of total weights compatible with the one obtained with 
the first one. However, the significance of the detection seems to be slightly 
lower. We are investigating whether this anomalous detection is due to some 
structures spatially localised or, on the contrary, the origin for such anomaly 
is homogeneously distributed in the sky. Our analysis indicates that the latter 
is the most plausible option. We are also analysing if the origin of the anomaly 
could be caused by a residual foreground contribution on the CMB map. If 
so, one should expect to find a frequency dependence for the anomaly, as the 
foregrounds themselves are frequency dependent. The same pattern of total 
weights has been found in the individual Q, V and W frequency channels of 
the first-year data. This might discard foregrounds as a possible origin of such 
anomaly. 



Other possible sources for explaining the anomaly are being considered. Firstly, 
the coincidence of the preferred direction (towards the ecliptic poles and 
the CMB dipole axes) naturally suggest possible unknown systematic effects 
fliil : liii ). Among them, errors in the beam pattern reconstruction are being 
studied. It is also interesting to note that the angular size of the mesh of the 
WMAP scan pattern defined by the combination of the spin and precession 
of the satellite is of the order of several degrees (js^. Second, the angular 
size at which the anomaly is detected is compatible with the size of primary 
CMB anisotropics due to topological defects such as texture fields (1511 1 or 
secondary anisotropics due to the Rees-Sciama effect ([s^. Alignment mech- 
anisms fi48i : i53i) were recently proposed which might be generalised to such 
structures. 
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